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The LHCb Collaboration has recently reported evidence for a CP asymmetry 
r-| approaching the percent level in the difference between D° — > tt + tt~ and D° — > 

K + K~ . We analyze this effect as if it is due to a penguin amplitude with the 
weak phase of the standard model c — > b — > u loop diagram, but with a CP- 
conserving enhancement as if due to the strong interactions. In such a case 
the magnitude and strong phase of this amplitude P^ are correlated in order 
to fit the observed CP asymmetry, and one may predict CP asymmetries for 
a number of other singly-Cabibbo-suppressed decays of charmed mesons to a 
pair of pseudoscalar mesons. Non-zero CP asymmetries are expected for D + — » 
K+K (the most promising channel for which a non-zero CP asymmetry has not 
yet been reported), as well as D° ->■ 7r°7r°, Df -> n + K°, and D s -> ir°K + . No 
CP asymmetry is predicted for D + — > 7r + 7r° or D° — > K°K° in this framework. 

PACS numbers: 13.25.Ft, 11.30.Er, 11.30.Hv, 14.40.Lb 



I Introduction 

Although CP violation was first observed in neutral kaon decays and CP asymmetries have 
been seen at the tens of percent level in B meson decays, the standard model describing 
these decays predicts naturally very small CP asymmetries in decays of charmed particles, 
of order 10~ 3 or less [U EJ |3]. These decays are dominated by physics of the first two 
quark families, with the contribution of the third family suppressed both by smallness of 
elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and by the relatively small 
b quark mass in the c — > b — > u penguin diagram. This is in contrast to the b — > t — > s 
penguin amplitude, which can profit from both a larger CKM factor and a much larger 
top quark mass. Following an early suggestion jl] that the penguin amplitude in D decays 
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may be enhanced by nonperturbative effects in analogy to the s — >■ d penguin amplitude in 
K — > 7T7T, recent studies [21 El [2] indicate that an order of magnitude enhancement is not 
impossible. 

The LHCb Collaboration has reported 3.5a evidence for CP-violating charm decays in 
the difference between CP asymmetries in D° — > K~ K + and D° — > 7r~7r + (6]: 

AA C P = A CP (K + K~) - A CP (n + n-) = [-0.82 ± 0.21(stat) ± 0.11(syst)]% . (1) 

For the decay of a charmed meson D to any final state / we are defining 

a m - r(z ^ /} ~ r( ^ /} (2) 

CPU) pp 4 J) + pp 4 J) • ^ 

Although the CDF II Collaboration at the Fermilab Tevatron does not see statistically com- 
pelling evidence for CP violation in either of these two decays, their results are consistent 
with those of LHCb [7]: 

A CP (D° K+K-) = (-0.24±0.22±0.09)%, A CP (D° tt+tt") = (0.22±0.24±0.11)% . 

(3) 

We calculate the corresponding 90% confidence level limits to be 
- 0.63% < A C p(D° K + K~) < 0.15% , - 0.21% < A CP (D° tt+tt-) < 0.65% . (4) 

The LHCb results have led to numerous hypotheses of physics beyond the standard 
model (e.g., [21 |8j [9j [10]) some of which had been studied earlier [11]. A more conservative 
approach, studying the above CP asymmetries within the Standard Model under relaxed 
assumptions about non-perturbative hadronic weak matrix elements, has been adopted 
recently in two papers applying flavor SU(3). Ref. [12] extended the hypothesis of triplet 
operator enhancement introduced in Ref. [I] by including in the effective weak Hamiltonian 
SU(3) breaking terms which are first order in the strange quark mass. A second work 
[T3] (appearing while we were writing up our results), applying a diagramatic flavor SU(3) 
approach similar to the one discussed by us below, associated H^-exchange and annihilation 
amplitudes with final state resonant effects [H]. While these two papers have some overlap 
with ours the specific assumptions and detailed predictions of the three studies, all based 
on flavor SU(3) analyses, are different. 

In the present paper we explore a scenario in which the standard model c — > b — > u 
penguin amplitude receives a sufficient enhancement from strong interaction physics to 
account for the effect observed by LHCb. This amplitude then must contribute to other 
direct CP asymmetries in decays of D°, D + , and to pairs of pseudoscalar mesons. Non- 
zero CP asymmetries are expected for D + — > K + T£° , as well as D° — > 7r°7r°, — > tt + K°, 
and D s — > tt°K + . No CP asymmetry is predicted for D + — > 7r + 7r° or D° — > K TC° in this 
model. The former receives no penguin contribution, being a AI = 3/2 process, while the 
latter involves a different penguin amplitude than the one we are considering. 

We perform the analysis in the context of a flavor-SU(3) model of charm decays pre- 
sented previously [151 [16]. We introduce notation in Sec. II and fit decay rates for singly- 
Cabibbo-suppressed (SCS) processes (including SU(3) breaking) in Sec. III. We then intro- 
duce a phenomenological penguin amplitude in Sec. IV to account for the CP violation 
observed by LHCb, and predict other CP asymmetries for SCS charmed meson decays. We 
summarize in Sec. V. 
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Figure 1: Relation between CF and SCS color-favored tree amplitudes. 



II Formalism and notation 

Cabibbo-favored (CF) charm decays may be characterized by amplitudes T, C, E, and 
A, corresponding to color-favored tree, color-suppressed tree, exchange, and annihilation 
flavor topologies [151 [16]. A fit to CF decays of D mesons to two pseudoscalar mesons leads 



to the following (\T\ > \C\) solution: 

T = 2.927 , (5) 

C = 2.337 e~ 1 15L66 ° = -2.057- 1.109 i , (6) 

E = 1.573 e il2a56 ° = -0.800 + 1.355 % , (7) 

A = 0.33 e i70A7 ° = 0.110 + 0.311 i , (8) 



quoted in units of 10~ 6 GeV. We note that complex conjugates of the amplitudes (E])-(JH1) 
give identical decay rates. To describe amplitudes corresponding to SCS processes the 
above amplitudes are multiplied by ±A, where A = tan^cabibbo = 0.2317. (Tiny phases of 
V u d, V cs , V c d, and V us are neglected. These contribute to negligible direct CP asymmetries 
in D° — > 7t + tt~ and D° — » K + K~ at a level below one per thousand by interference of the 
above amplitudes with SU(3)-breaking penguin amplitudes discussed in Section III.) The 
relation between CF and SCS color-favored tree amplitudes is illustrated in Fig. [TJ 

In order to account for SU(3) breaking in the SCS T amplitude we may use the following 
expressions: 

Td°-^-k+-k- = Td+^tt+tt = T D +^ w+K o = T n , (9) 
T d o^k+k- = T D +^ K +K° = T K , (10) 
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Figure 2: Relation between CF and SCS A amplitudes. 



where, neglecting the contribution of f-(q 2 ) at q 2 = m 2 K , 

n If / 9 \ I " 9 9 > \ / 

\f + (D°->K-){mi)\ m z D - m l K 

^ _ rp \f+(D ^K-)(m 2 K)\ Ik / 1n x 
\f + (D°->K-){mi)\ U 

Similarly, the relation between CF and SCS A amplitudes is illustrated in Fig. [2j Here, 
one may introduce SU(3) breaking as follows: 

Jd+ 

A d +^ +k0 ~ ^ D +^ +7r o = A . (14) 

We know the relevant decay constants [17] and meson masses [18] (in GeV): 

U = 0.13041; f K = 0.1561; / D+ = 0.2067; f D + = 0.2575; (15) 
m D o = 1.8648; m w = 0.13957018; m K = 0.493677. (16) 

The following approximate values are also known for the form factors from semileptonic 
D° decays PHED]: 

|/+(do^-)K)|~ 0.705, (17) 
\f +{D ^ K -){ml)\^ 0.768, (18) 
\f+^-,K-){m 2 K )\^ 0.811. (19) 

After using relevant form factors and decay constants, we find, in units of 10~ 6 GeV, 

T w = 2.87 , (20) 
T K = 3.70, (21) 
A D+ = (0.89 + 2.50 i) x 10" 1 . (22) 
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Figure 3: Penguin diagrams leading to a non-zero amplitude P = Pd + P s in the presence 
of imperfect cancellation between intermediate d and s quarks. 




Figure 4: Penguin annihilation diagrams leading to a non-zero amplitude PA = PAd + PA s 
in the presence of imperfect cancellation between intermediate d and s quarks. 

The tree-level amplitudes for D -» 7T+7T-, D° K + K~, and D° 7r°7r° involve the 
following respective combinations (in units of 1CT 7 GeV): 

-X(T V + E) = -4.80- 3.14 i, (23) 
\(T K + E) = 6.72 + 3.14 i, (24) 
\(C-E) = -2.91- 5.71 1, (25) 

as well as SU(3)-breaking terms which we shall now introduce. 

Ill Fits to decay rates including SU(3) violation 

The penguin amplitude P for c — > u transitions is normally thought to be very small because 
the contributions of d and s quarks in the intermediate state cancel one another [21J. If 
we regard this cancellation as inexact due to SU(3)-violating masses of intermediate-state 
particles, we can regard the penguin amplitude P as a proxy for SU(3) violation (see Fig. 
[3]). It will then have the same weak phase as other standard model contributions to D 
decays. The same can be said for a penguin annihilation (PA) amplitude, contributing 
only to D° decays. It corresponds to the exchange processes cu —> ss and cu —> dd followed 
by ss or dd annihilation into a pair of charge-conjugate pseudoscalar mesons (Fig. H]). 

The amplitudes for D° — > 7r + 7r~, D° — > K + K~~, and D° — > 7r°7r° then may be expressed 
as shown in the first three lines of Table [B Given the magnitudes of the relevant amplitudes 
determined from the decay rates [T6], [22], in units of 10 -7 GeV, 

\A(D° -> = 4.70 ±0.08, (26) 
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Table I: Representations and comparison of experimental and fit amplitudes for SCS decays 
of charmed mesons to two pseudoscalar mesons. 



Decay 


Amplitude 


\A\ (10- 

1 IV 


1 GeV) 


x 2 


Mode 


representation 


Experiment 


Theory Fit 




D° -> tt+tt- 


-A (T ff ± P) ± (P ± PA) 


4.70±0.08 


4.70 





D° K + K~ 


A (T K ± P) ± (P ± PA) 


8.49±0.10 


8.48 


0.01 


D° -> ttV 


-A (C - P)/V2 - (P + PA)/y/2 


3.51±0.11 


3.51 





P> + 7T + 7T 


-\(T W + C)/V2 


2.66±0.07 


2.26 


33 


D° -> 


-(P + PA) + P 


2.39±0.14 


2.37 


0.02 


P>+ -> iT+K 


X(T K -A D+ ) + P 


6.55±0.12 


6.87 


7 


P> + -> TT + PJ 


-A (T^ - A) + P 


5.94±0.32 


7.96 


40 


P>+ ->• 7T°ir + 


-A(C + A)/V2-P/V2 


2.94±0.55 


4.44 
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\A{D°^K + K-)\ = 8.49 ±0.10, (27) 
V2\A(D° 7r°7r°)| = 4.96 ±0.16, (28) 

one may then plot circles with these radii and centers defined by Eqs. ( I23H251) to solve 
for a common value of P + PA. The existence of a self-consistent solution for P ± PA is 
supported by a ^-minimization fit, which leads to 

P+PA = [(0.44±0.23)±(1.41±0.36) i] x KT 7 GeV ; X 2 /d-o.f. = 0.012/1 = 0.012 . (29) 

The construction and the corresponding A% 2 = 2.3 error ellipse (corresponding to 68% 
probability) are shown in Fig. [5j 

Using the extracted value of P+PA we apply a similar construction technique to extract 
P. The relevant parts of the tree-level amplitudes that determine centers of the circles are 
as follows (in units of 10 GeV): 

- (P + PA) = -0.44 - 1.41 i , (30) 

\(T K -A D+ ) = 8.37- 0.58 i, (31) 

-A (T w — A) = -6.40 + 0.72 i , (32) 

\(C + A) = -4.51- 1.85 i. (33) 

The relevant experimental rates [Ml lead to amplitudes once again determining the 
radii of the circles as follows (in units of 10 -7 GeV): 



\A(D°^K°TC)\ = 2.39 ±0.14, (34) 

\A(D + ->■ K + K°)\ = 6.55 ±0.12, (35) 

\A(D+ -> P°7r+)| = 5.94 ±0.32, (36) 

V2\A(D+ -»■ K + 7t°)\ = 2.94 ±0.55. (37) 

X 2 -minimization gives us 

P= [(-1.52 ±0.15) ±(0.08+°;|) i] x 1(T 7 GeV ; X 2 /d.o.f. = 54/2 = 27 . (38) 
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-10 -5 5 10 -0.8 -0.6 -0.4 -0.8 

Re(P + PA) (10~ 7 GeV) -Re(P + PA) (1Q~ 7 GeV) 

Figure 5: Construction to determine P + PA. The relative sign between the left-hand and 
(magnified) right-hand panels is due to the fact that the vector P + PA points toward the 
origin in the left-hand figure. 




-10 -5 5 10 1.0 1.2 1.4 1.6 1.8 3.0 



Re(P) (10 7 GeV) -Re(P) (10 7 GeV) 

Figure 6: Construction to determine P. The relative sign between the left-hand and 
(magnified) right-hand panels is due to the fact that the vector P points toward the origin 
in the left-hand figure. 
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The construction and the corresponding 68% error ellipse (Ax 2 = 2.3) are shown in Fig. [6j 
In Table [I] we quote the representations and compare the experimental and fit amplitudes. 
Using the extracted values of P and P + PA we find 

PA = [(1.95 ± 0.38) + (1.34 ± 0.71) i] x 10~ 7 GeV . (39) 

We recall the two-fold ambiguity permitting amplitudes which are complex conjugates of 
d30J)-P3D and (EHD-([39D. 

The poor \ 2 i n this fit is driven primarily by the large contribution from the — » tt + K° 
amplitude. It is quite possible that our description of SU(3) breaking in this quantity is 
imperfect. In any case, the large experimental errors on the SCS decays of D s to two 
pseudoscalar mesons will hinder the study of CP-violating asymmetries in their decays for 
some time to come, so we shall not be greatly concerned with such decays for the present. 

IV Description and prediction of observed direct CP 
asymmetries 

We now consider the effects of an additional phenomenological-penguin amplitude P b , the 
weak phase of which differs from the weak phase of P and PA by the CKM-angle 7. (The 
subscript b refers to a b quark in the intermediate quark loop in Fig. |3j) In Table [III we 
summarize the amplitudes for SCS processes obtained in the previous section, and extend 
the amplitude representations to include P b . The quantities 0^ = ArgfTj] denote the strong 
phases of the non-P fe amplitudes with respect to T. (The amplitudes Tf include factors 
±A.) 

Table II: Fit amplitudes for SCS charmed meson decays including P and PA, and their 
representations including P b . 



Decay 


Amplitude 


4 = Aig[T f ] 


mode 


representation 


(degrees) 


D° 7T+7T 

D° -»■ K+K- 
D° -> 7r°7r° 


-A (Tjr + E) + (P + PA) + P b 
X (T K + E) + (P + PA) + P b 
-A {C - E)/y/2 - (P + PA) - P b /V2 


-158.5 
32.5 
60.0 


D+ 7T+7T 


-X(T^ + C)/V2 


126.3 


D° ->■ K°K" 


-(P + PA) + P 


-145.6 


D + -> K + K° 


X(T K -A D+ ) + P + P b 


-4.2 


D+ -> n+K° 


-A (T n - A) + P + P b 


174.3 


D+ -± TT°K + 


-X (C + A)/V2 - P/V2 - P b /V2 


16.4 



In general, the amplitude for D — > f may be written as follows: 

A(D -> /) = \T f \ e 1 ^ (l + r /e i(7+</,/) ) , (40) 

where T$ represents terms that have the same weak phase as the tree-level terms contribut- 
ing to that amplitude, 0^ represents the strong phase of Tf, Tf represents the ratio of the 
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magnitude of the CP- violating penguin contribution to that of Tf, 7 represents the weak 
phase of the CP-violating penguin (it is the same as the CKM angle), and <f>* is the strong 
phase of the CP- violating penguin relative to Tf. Let us take the example of the process 
D° — > 7r + 7r~ for clarity. Then 



67T^7r 
T 



T W+7T - = -\(T W + E) + (P + PA) , (41) 

= ArgfT^-] , (42) 

- J^L, (43) 

|-^7r+7T I 

= Arg[P b ] 7 . (44) 



The amplitude for D — > / may be written as follows: 

A(D J) = [Tfle^r (l + r/e*^^) . (45) 

For a two-body decay, the rate is proportional to the absolute square of the amplitude. 
Thus, one may now define a CP asymmetry as follows: 



Acp(f) 



T(D^f)-T(D^f) 

r(£> —>■/) + r(£> — >■ /) 

2 r ^ sin 7 sin <\J 
1 + r j + 2 r/ cos 7 cos <^ 

2p|T/-| sin 7 sin(5 — 
\T f \ 2 + p 2 + 2 p \Tf\ COS7 cos(<5 — 4>t] 



(46) 



where in the final step we have used P& = pe* ( - <5+7 * ) . 

The LHCb result ([T]) [6] may be used as a constraint on the magnitude and strong phase 
of the CP-violating penguin P^. [Note added in proof: To lowest order in p, all asymmetries 
Acp{f) depend on the combination psin7. Thus, if we impose the AAcp constraint, our 
predictions for other asymmetries are the same for any weak phase of Pb as long as effects 
of higher order in p are negligible. We thank N. Deshpande for a question leading to this 
result.] We use the following relationships: 



\Tk+k-\ +P + 2 P \Tk+k-\ cos7Cos(<) - (pf ) 

a 2p\T w+n ~\ sin7sin(^-^ + ^) 

A CP {n tf ) = -— — ■ (48) 



We may use the theory fit results quoted in Tabled] for \Tf\. The strong phase can be 
taken from the results quoted in Table HT1 The CKM angle 7 may be taken to be 77° [18J. 
Corresponding to each value of S allowed by the AAcp constraint, one may extract the 
allowed values of p. In addition we expect |P&| < \Tf\ (r/ < 1), which in turn restricts us 
to small values of p. In Fig. [7] we plot the allowed values of p as a function of 5 using Eqs. 
(14T|) and (14"8"|) . For a wide range of 5, a penguin amplitude of magnitude 0.01 x 10~ 7 GeV, 
or (9(0.1%) of the D° — > K + K~ amplitude, is sufficient to account for the observed value 
of AAcp- This is in accord with a conclusion reached in Ref. [3]. 
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Figure 7: p and 5 allowed by the measured range of AAqp- The (red) line represents 
the central value, while inner (blue) and outer (green) bands respectively represent 68% 
confidence level (lcr) and 90% confidence level (1.64a) regions based on error in AAqp- 
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Figure 8: Acp(K + TC°) as a function of the allowed values of 5. The (red) line represents 
the central value, while inner (blue) and outer (green) bands respectively represent 68% 
confidence level (la) and 90% confidence level (1.64a) regions based on error in AAqp- 
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Figure 9: Acp as a function of the allowed values of 5. The (red) line represents the central 
value, while inner (blue) and outer (green) bands respectively represent 68% confidence 
level (la) and 90% confidence level (1.64cr) regions based on error in AAcp. 
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The constraint on p as a function of 5 now allows us to predict CP asymmetries in other 
channels such as D + — > K + K as a function of the angle 5. In Fig. [8] we plot Acp(K + K ) 
as a function of 5. Values of p are plotted only for the range of 5 consistent with the limits 
(jll), which will be specified shortly. In Fig. [9]we plot Aqp for the final states K + K~ , 7r + 7r~, 
and 7r°7r°. The limits (jl]) imply the following allowed ranges of 5: 

A C p(K + K-) => 0.50 < S < 3.57 , A CP (7i + n~) 0.49 < S < 3.57. (49) 

Figs. [7J and [8] are plotted only for values of 5 consistent with both these limits. Note the 
correlation between the CP asymmetries in the channels D° — > tt tc and D + — > K + K° . 
More precise measurements of the individual asymmetries in D° — > 7r + 7r~ and D° — > K + K~ 
can help to pin down the unknown strong phase S. 

As mentioned in the preceding section, all the contributions to Tf listed in Table HT1 
involve an ambiguity due to complex conjugation. Thus, the phase $ T has a sign ambiguity, 
(f)j, — >• — 4>j>, which is common to all final states /. The CP asymmetry (14"6"|) is approximately 
invariant under a joint transformation, 

4> S T -» -<p T , 5 -> 7r - 5 , (50) 

neglecting a very small contribution to the asymmetry quadratic in p/\Tf\. Thus, while 
plots similar to Figs. [71 [8] and [9] may be plotted with 5 — > n — 5, the correlations between 
asymmetries in different decay modes are invariant under this redefinition. 

We have left out decay asymmetries since the corresponding branching ratios have 
large fractional errors. The process D + — > 7r + 7r° does not depend on Pj, in the isospin 
symmetry limit, and therefore its CP asymmetry is zero at this high level of approximation. 
The CP asymmetry in D° — > K°K depends only on a penguin annihilation diagram as 
there are no u quarks in the final state. If it is found to be non-zero, our discussion must 
be expanded to include the possibility of CP violation due to interference between a (PA)b 
amplitude involving a b quark in the loop and an SU(3) breaking term in E. 

V Discussion and summary 

The observation by the LHCb Collaboration of a difference between the CP-violating asym- 
metries in D° — > K + K~ and D° — > 7r + 7r - likely implies observable asymmetries in other 
decays of charmed mesons to pairs of pseudoscalar mesons. The present description of 
that difference assumes that a penguin amplitude with an intermediate b quark, normally 
thought to provide a contribution below the observed effect, is amplified by CP-conserving 
physics (e.g., unforeseen QCD effects) to an extent which can account for the asymmetry. 
In that case several direct CP asymmetries are predicted as functions of a single strong 
phase difference 5. These include asymmetries in the individual decays D° — > K + K~ and 
D° — > 7r + 7r~, as well as D° — > 7r°7r° and D + — > K + K . These asymmetries are typically 
of order (a few) xlO -3 , and the latter two are correlated with one another. Experimental 
limits (jl]) on the direct CP asymmetries in D° — > K + K~ and D° — > 7t + tt~ [7\ provide 
constraints on 5. The observed asymmetry [23] Aqp(D + — > K + K°) = (7.1 ± 6.1 ± 1.2)% 
carries far too large an uncertainty at present to test its prediction. [Note added in proof: 
(1) The CDF Collaboration has now reported a value of AA CP = (-0.62 ± 0.21 ± 0.10)% 
|24j. (2) We thank Anze Zupanc for reminding us that the Belle Collaboration has reported 
the much more precise value A CP (D + K + K°) = (-0.16 ± 0.58 ± 0.25)% [25].] 
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In Fig. El while Acp(K + K~) and Acp{it + ii~) are predicted to have opposite signs 
for a wide range of 5, their relative magnitudes provide information about 5, with the 
ratio \Acp{it + ii~) / Acp(K + K~)\ exceeding 1 for the mid-range of 5 and behaving as a 
decreasing function of 5. Thus, better measurements of these individual asymmetries will 
enable improved predictions of asymmetries such as Acp(K + K ) and Acp(7r°7r )- We look 
forward to improvement of many of these determinations. 
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